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ABSTRACT 



We obtain the general forms for the current density and the vorticity from the integra- 
bility conditions of the basic equations which govern the stationary states of axisymmetric 
magnetized self-gravitating barotropic objects with meridional flows under the ideal MHD 
approximation. As seen from the stationary condition equations for such bodies, the presence 
of the meridional flows and that of the poloidal magnetic fields act oppositely on the internal 
structures. 

The different actions of these two physical quantities, the meridional flows and the 
poloidal magnetic fields, could be clearly seen through stationary structures of the toroidal 
gaseous configurations around central point masses in the framework of Newtonian gravity 
because the effects of the two physical quantities can be seen in an amplified way for toroidal 
systems compared to those for spheroidal stars. The meridional flows make the structures 
more compact, i.e. the widths of toroids thinner, while the poloidal magnetic fields are apt to 
elongate the density contours in a certain direction depending on the situation. Therefore the 
simultaneous presence of the internal flows and the magnetic fields would work as if there 
were no such different actions within and around the stationary gaseous objects such as ax- 
isymmetric magnetized toroids with internal motions around central compact objects under 
the ideal MHD approximation, although these two quantities might exist in real systems. 
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1 INTRODUCTORY ANALYSIS AND MOTIVATION 

1.1 Theoretical treatment of stationary states of 

axisymmetric magnetized self-gravitating barotropes 
under the Ideal MHD approximation 

In this paper, we give new expressions for the current density 
and the vorticity vector inside the stationary and axisymmetric 
magnetized self-gravitating barotropes with internal gaseous mo- 
tions under the ideal magnetohydrodynamics (MHD) approxima- 
tion. Although stationary states of the magnetized self-gravitating 
barotropes without internal flows have been investigated in many 
old paper s (e.g., 'Chandrasekhar & Fermi 1953; Liist & Schliiter 
19541; [Feri -aro 1954; GieOestad 1954; Roberts 1955; Chandrasekhar 
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pers (e.g.,|T omimura & Eriguchi 2005; Yoshida & Eriguchij l2006l ; 
lYoshidaetalll200a; lOtani et alj 120091 ; IFuiisawa et all 120 la) , the 
effects of the internal flows on structures of magnetized self- 
gravitating barotropes have barely studied so far. 

Dynamic equilibrium equations for stationary states of mag- 
netized self-gravitating bodies with internal flows are given by 



-Vp = -V<;/>g - V0C 



-V|v|- 



D X oj H j X B . 

cp 



(1) 
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where p,p, 4>g, 4>c,v,(jJ,c, j, and B are the density, the pressure, 
the gravitational potential of the gaseous body, the gravitational po- 
tential of external objects, the fluid velocity, the vorticity vector, the 
speed of light, the current density, and the magnetic field, respec- 
tively. In this study, the gravitational potential of external objects is 
assumed to be given by 

<Pc = , (2) 
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where Mc, G, and r denote the mass of the central external object, 
the gravitational constant, and the distance from the central object, 
respectively. Using the expressions for the current density and the 
vorticity vector inside the stationary, axisymmetric and infinitely 
conducting barotropes (for details, see AppendixjA), we may sim- 
plify the forth term in the right-hand side of Eq. (|2^as 

f 1 dQ('i) 



V X uj — • 



A-Kcp d\I/ 

»'(Q)V(5, 



(forB = 0, Tip ^ 0) 
(foiB = 0,i>p = 0) 



where ^ and Q are the flux and stream functions, respectively (see 
Eqs. JAU and iA2t '). Q, and v are arbitrary functions of the given 
argument, and R and e.^ are the radial coordinate and base vector 
for the cylindrical coordinate with respect to the symmetry axis, 
respectively. We describe the details of these arbitrary functions 
in Appendix |A] and Appendix |B] Here, Vp and v^ stand for the 
poloidal velocity and the toroidal velocity, given by 



4^^«^^^' 



= RQ.+ 



1 rfQ(^) 
47rcp d* 



B^ 



(4) 



(5) 



with tp and B^ being the azimuthal angle and the toroidal mag- 
netic field, respectively. The stream function Q is a function of the 
flux function ^ if and only if B 7^ and Vp 7^ 0. Thus, equi- 
stream function surfaces are similar to equi-magnetic flux function 
surfaces if B 7^ and Vp 7^ 0. Similarly, the Lorentz force term 
appearing in the fifth term in the right-hand side of Eq. ([T) may be 
rewritten as 



-J X B = 



< 47r 

I 0, 



1 dQ(*) 



cp d* 



X B , (for B 7^ 0) 
(for B = 0) 



such effects have not been investigated because stationary states of 
axisymmetric magnetized barotropes with meridional flows have 
not been obtained. Thus it is one of the puiposes in this study to 
show the above statements could be proved to hold in some numer- 
ical examples. 

(3) 1.2 Toroids: Best astrophysical systems in which both 
meridional flows and magnetic fields would work 
simultaneously but differently 

In order to show clearly the oppositely working effects of the above 
mentioned two quantities, meridional flows and poloidal magnetic 
fields, we will show numerical results of stationary configurations 
of axisymmetric magnetized self-gravitating toroidal barotropes 
with meridional flows which locate around central point masses in 
the framework of Newtonian gravity under the ideal MHD approx- 
imation. For toroidal configurations, matter distribution changes 
in wider space would be expected compared to size changes of 
spheroidal objects because toroids could change their shapes to two 
opposite directions, i.e. to the outside direction and to the inside di- 
rection of the toroids. Thus we will solve stationary states of ax- 
isymmetric magnetized barotropic toroids with meridional flows 
under the ideal MHD approximation and clarify the oppositely 
working effects explicitly in this study. 

Concerning self-gravitating toroidal configurations or disks, 
we need to take recent results of fully general relativis- 
tic (GR) numerical simulations into account. These simula- 
tions show that a few-solar-mass black hole and a highly 
dense toroid whose maximum density can reach 10^*^ - 
IQ^^g/cm around the blac k hole could be form e d after merg 



dp 
P 



where At(*l/) is an arbitrary function. From Eqs. (1T]-(|6], we obtain 
Bernoulli's equations for the present situation: 

GMc 1 „ , /--I' 

-<^g + {\llp\^ + l.^) + flt,^n(*) + / A.{*)d* + c , 

r 2 J 

(for B 7^ 0) 

GM^ 1 , , /-Q 

-<*g + — -(t^p|2+„2)+y i,^Q)dQ + C, 

(forB = 0,i>p # 0) 

GA/c /■« 2 

-4>g + h/ BndK+C, 

(forB = 0,i>p = 0) 

where C is an integration constant. It should be emphasized that 
these three expressions cannot be united to a single expression as 
can be understood from Eqs. (|3) and (|6). Note that in this study, 
we do not consider the case of Bp — 
Bernoulli's equations differ from Eq. (|7). 

By comparing Bernoulli's equation for the case of i? = 
and Vp = with that for the case of B = and Vp 7^ 0, we 
may expect that the presence of the meridional flows tends to de- 
crease the volume occupied by the fluid. This is because the ki- 
netic energy of the meridional flow contributes to the fluid as posi- 
tive rampressures (dynamic pressure), which result in reducing the 
fluid pressure and consequently decreasing the density of station- 
ary fluid objects. In contrast to the effect of the meridional flow, for 
magnetized equilibrium configurations with Dp = whose mag- 
netic fields are generated by positive toroidal currents, the poloidal 
magnetic field is apt to expand the fluid region to the direction per- 
pendicular to the symmetry axis like the centrifugal force (see, e.g. 
old papers Ferraro 1954; Gjellestad 1954 ; Roberts 1955; and re- 
cent paper s Tomimura & Eriguchi 2005; Yoshida & Eriguchil200q ; 
lYoshida et al.i2006 ; iFuiisawa et al.i2 012). Therefore, simultaneous 
presence of the meridional flows and the poloidal magnetic fields 
could result in almost no changes in the matter distributions. So far. 



ing o f binary neutron star s (IShibata & Urv u 2000; Shiba ta et al 
20031 ; IShibata et al] |2005| ; iKiuchi et alfHoOa ; .Hotokezaka et al 
201 ih . after mergi ng of a neutron star and a black hole 



in binary s ystems ('Shibata & Urvii' '2006':'Shibata & Urvu' '2007| 
Shibata & T aniguchi 2008; Kyutoku et al. 2010; Kyuto ku et al 
2011b or after collapsing of a supermassive rotating star (Shibati 



woo"; 'Shibata & ShapirJ2 002l;IShibatal 2003'; 'Sekiguchi & Shibad 
2004; Sekiguchi & Shibatj bOOTT Tsekiguchi & Shibata 201^ 
Therefore, dense toroids and central compact objects could be 
formed after collapsing or merging of compact objects. Similar 
kinds of systems with magnetic fields have a lso been investigated 
by se v eral groups (e.g . 'Naravan et al.l uOOlt IShibata & Sekiguchil 
I2OO5I ; lDuezetal.l lm)6; Shibata et all |2007|) . Although, in or- 
der to understand the origin and dynamical formation pro- 
cesses of these systems, we must take into account many re- 
alistic physics and compute stationary configurations with mag- 
netic fields in GR, nobody has yet succeeded in solving sta- 
tionary states both with poloidal and toroidal magnetic fields 
in GR at present. Therefore, we explore such stationary states 
of axisymmetric magnetized barotropic systems in the frame- 
work of Newtonian gravity. Although there were many papers 
to obtain magnetized stationary states of disks/toroids only with 
poloidal magnetic fields (e.g. Bisnovatyi-Kogan & Blinni kov 197j; 



Bisnovatvi-Kogan & Seidffsj 1 19851 ; iBaureis et all 1 19891 ; iLi & Shi] 



199g) and disks/toroids only with toroida l magnetic fields in- 



side the disks/toroids (e.g. Okada et al.lll989l ; lBaneriee et al.lll995l ; 



iGhanbari & Abbassil 120041) . no solutions both with poloidal and 
toroidal components of magnetic fields have been obtained yet. 
This is because it has been difficult to solve the Grad-Shafranov 
equation as well as the equations of motion consistently by some 
means. Concerning this type of systems, the most general formula- 
tion was derived bv lLovelace et al.l ( 119860 systematically. However, 
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iLovelace et alj ( 119860 computed solutions only with poloidal mag- 
netic fields for non-self-gravitatin g disks. 

Recently, Otani et al. ( 2009) have obtained magnetized self- 
gravitating equilibrium states both with poloidal and toroidal mag- 
netic fields self-consistently in t he framework of Newt o nian g rav- 
ity. Their method is based on iTomimura & Eriguchil ( 2005h. In 
this st udy we have extended the method employed in lOtani et alJ 
( 12009 ) to the most general configurations for the stationary states 
of axisymmetric magnetized barotropic toroids with meridional 
flows under the ideal MHD approximation and obtained se- 
quences of stationary states. Comparing these results w ith those 
of non-magnetiz ed toroids without meridional flows (e.g. lOstriken 
ll964llWonal974|), with those of magnetized toroids without merid- 
ional flows (e.g. lOtani et alj2009.) or with those of non- magnetized 



toroids with meridional flows (e.g. lEriguchi et al.ll986l) , we will be 
able to clearly see the effect of the presence of both physical quan- 
tities, i.e. meridional flows and magnetic fields as explained in the 
previous subsection. 



2 BRIEF DESCRIPTION OF THE PROBLEM 

In this study, as mentioned, we investigate stationary configura- 
tions of axisymmetric magnetized polytropic toroids with internal 
fluid motions. We consider inviscid and infinitely conductive toroids 
with equatorial symmetry located around central point masses in 
the framework of Newtonian gravity. Since a si milar problem, bu t 
without meridional flows, was already treated bv lOtani et al.l ( 120091) 
and our strategy is basically the same as theirs, we briefly summa- 
rize the basic equations, boundary conditions, and solving scheme. 
We describe the details of physical quantities and dimensionless 
forms in Appendix Icl and the numerical scheme in Appendix IdI 

Continuity and pressure balance equations for stationary states 
are respectively given by 



V • {pv) = , 



(8) 



1^ ^, ^fG^c 

-Vp = -V(/.9 + V 

P \ r J 2 



V|i;^l + vxuj + -(-xB 



The Poisson equation for (pg is 

A(t)g = 4TTGp. (10) 

We make use of the polytropic equation of state 



P = Kop 



(11) 



where Ko and A^ are a constant and the polytropic index, respec- 
tively. 

Maxwell equations for stationary electromagnetic fields are 



V B = 0, 

V X E = 0, 

V X B = 47r^ 



(12) 
(13) 
(14) 



where E is the electric field, which is determined by the perfect 
conductivity condition 



E+ - X B ^0. 
c 



(15) 



The magnetic flux function '3/ is, in terms of the azimuthal compo- 
nent of the vector potential A^p, defined by 



From the azimuthal component of the Maxwell equation ( I14t . we 
obtain 



R 



dR 



Rim 



+ 






'AtvR 



3-fi 



(17) 



where the cylindrical coordinates (R, ip, z) have been used. This 
equation is equivalent to the so called Grad-Shafranov equation for 
this problem, but we treat the right-hand side of Eq. ( 117) as the 
source term of the differential operator in the left-hand side of Eq. 
( 1171 ) even though j^ includes the term proportional to the left-hand 
side of Eq. J17| ( if Vp 7^ 0. Explicit expression for j^ is given in Eq. 
([at}. In order to solve Eq. ( I17t numerically, it is useful to transform 
Eq. Ml\ into 



A(yl^ sin ip) = — 47r— sin Lp 



(18) 



where A denotes the Lap lacian (see, e.g., iTomimura & Eriguchil 
I2OO5I and lOtani et alll2009l) . 

Because isolated mass and current densities are considered in 
this study, boundary conditions at infinity for the gravitational po- 
tential (j)g and the vector potential A^ are respectively given by 



O 



A^ 



O 



(19) 



Solutions for the present problem are obtained by solving the pres- 
sure balance equation and the two Poisson equations for (jig and 
A^ sin ((3 with the boundary condition ( I19t . In order to impose 
the boundary condition ( 119) automatically, in this study, the two 
Poisson equations dlOt and l ll8t are converted into the two integral 
equations, given by 



= -G 



A^(r,G)wi^p 



-p-^ — —smip d r 
c I \r — r'\ 



(20) 



(21) 



^ = RA^ 



(16) 



Integrating Eq. (|9), we obtain Eq. (|7}. Eqs. (|7}, ( I20t , and ( I21t are 
the master equations for the present proble m, which are s olved with 
a variant of the numerical scheme used by lOtani et al.l ( |2009|) (for 
(9) details, see Appendix iDt. 

In general, the toroid cannot approach indefinitely to the cen- 
tral object because gravitational effects of the central object can 
unboundedly increase as the distance from the central object to the 
toroid decreases and any forces counteracting the gravity cannot 
stanch the matter flow shedding from the inner edge of the toroid 
if their distance is shorter than some critical value. In the present 
numerical scheme to obtain magnetized toroids, the distance from 
the central object to the inner edge of the toroid (or the width of the 
toroid) is characterized by a dimensionless parameter q, defined by 

'7=-^ , (22) 

Pouter 

where -Rinnor and -Router are the shortest and the longest distances 
from the symmetry axis to the toroid, respectively. In teiTns of q, 
this disappearance property of the equilibrium sates describes as 
the existence of qc such that there is no stationary solution of the 
toroid for q < qc. Note that the value of q^ depends on what param- 
eters characterizing equilibriu m sequences keep c onstant when the 
value of q changes. Following lotaiii et al.l ( l2009f) . we call equilib- 
rium solutions characterized by q = qc the critical configuration or 
the critical state. The distance from the symmetry axis to the inner 
edge of the toroid for the critical configuration is named the critical 
dist ance. In this s tudy, we focus only on the critical configurations. 
lOtani et al.l ( 120090 investigated the critical configuration for 
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the magnetized toroids without meridional flows and found the fol- 
lowing properties, (i) The critical configuration features cusp-like 
structures at the inner edge of the toroids. (ii) The critical config- 
uration rotates very slowly. This implies that the critical toroids 
are mainly sustained by the balance among the magnetic forces, 
the gravity of the central objects and the pressure gradients, (iii) 
The critical distances are almost independent of the mass ratio of 
the toroids to the central objects, (iv) The critical distances are 
much larger than QGMc/c? , the radius of the innermost stable 
circular orbit for the Schwarzschild black hole with the gravita- 
tional mass Mc- This means that making use of Newtonian gravity 
is reasonab l e to in vestigate structures of the toroids considered in 
lOtani et alj i2009l) . 



3 NUMERICAL RESULTS 



Following lOtani et al.l (2009), we consider two polytropic indices 
TV = 1.5 and N = 3 only in the present study. As for the ar- 
bitrary functions, which need to be specified to obtain particular 
solutions, we employ the same func t ional forms as those used in 
lYoshida et al] I2OO6) and Otani et al.l J2009l) except for the toroidal 
current function ji. For the toroidal curre nt function fi, we cho ose 
the same functional form as that used in iFuiisawa et al.l ( 120 12h in 
order for the inner edge of the toroid to have stronger magnetic 
fields. As for the stream function Q\ which does not appear in 
lOtani et alj i2009l) . a similar functional form to that of the poloidal 
current function k is employed. Details of the functional forms of 
the arbitrary functions are collected in Appendix IA2I 

As mentioned before, recent numerical simulations performed 
with numerical relativity show that geometrically thick toroids ro- 
tating around black holes form after mergers of neutron star bi- 
naries, mergers of black hole-neutro n star binaries, or collapses 
of supermassive rotating stars (e.g. ISekiguchi & Shibatal uOOTI : 
ISekiguchi & Shibatal201lh . Typical values of physical quantities of 
such black hole-toroid systems are given by Mt — 1.0 x 10^^ A4q, 
Mc = 5.OM0 (Mt/Mc = 2.0 X 10"^), and pc = 1.0 x 
lO^^g/cm^, where Mt and pc are the mass and the maximum 
density of the toroid, respectively. Since these black hole-toroid 
systems are of large significance in high energy astrophysics, we 
focus on the models characterized by the mass ratio Mt/Mc = 
2.0 x 10""^ and use these values of the mass of the central object 
and the maximum density of the toroid to estimate values of other 
physical quantities with physical dimension, e.g. \B\ and Rinner- 

To check numerical accuracies of the stationary configurations 
obtained in this study, we estimate values of a virial relation, which 
vanishes for exact stationary solutions: 



VC = 



2T + W + 3U + m 



W 



(23) 



where T,W,Il, and OT are the kinetic energy, the gravitational en- 
ergy, the volume integral of the pressure and the magnetic energy, 
respectively (for details, see Appendix |C}. As shown later, all the 
stationary configurations obtained in this study are within accept- 
able accuracy (for details, see Appendix Et. 



3.1 Widening of the widths of toroids: Effect of the localized 
poloidal magnetic fields 

The counter effects of the meridional flows against the magnetic 
forces on structures of magnetized toroids would be clearly seen for 
toroids with rather widened shapes due to poloidal magnetic fields. 



Thus, in this subsection, we will try to compute magnetized toroids 
with highly localized poloidal magnetic fields because such equi- 
librium configurations could be toroids with a rather small value of 
q, i.e. th e width of toroids on the equatorial plane being rather wide 
(see e.g. lFuiisawa et alj2012h . 

To investigate effects of the localized poloidal magnetic fields 
on the toroid structures, no fluid flow inside the toroid is considered 
here. Thus, values of the parameters Qo, ^q, ct, and d are taken to 
be Qo — 0, CIq = 0, a — 0, and d = (for details, see Appendix 
lA2l (. As for Ko, following Fotani et all J2009l) , we take kq = 4.5. 
Since there are no rotation and no meridional flows, the toroids are 
in stationary states by the balance among the gravitational force of 
the central object, the Lorentz force and the gas pressure gradient. If 
there could be very strong poloidal magnetic fields near the central 
objects, stronger gravitational forces of the central objects could be 
balanced by the strong magnetic forces near the central objects. For 
such situations the toroids could be elongated toward the central 
obj e cts and have wider wi dths . 

iFuiisawa et al.l ( 120121) showed that the poloidal magnetic field 
distributions substantially depend on the parameter m in the arbi- 
trary function fJ,{^) for magnetized stars. In particular, they found 
that negative values of m result in concentration of the poloidal 
magnetic fields near the symmetry axis of magnetized stars. Thus, 
it is expected that we may obtain toroids in which poloidal mag- 
netic fields are concentrated near the inner edge of the toroids by 
choosing an appropriate value of m. 

We show our numerical results for critical configurations of 
A*' = 1.5 and A'^ = 3 polytropic sequences in Table [T] and density 
contours on the meridional cross section of A'^ = 3 polytropes with 
m = —1.4, m = —0.7, m = —0.1 and m = 0.5 in Fig.[T] As 
seen from these table and figure, the critical distance decreases as 
the value of 771 decreases. The density distribution of the m = —1.4 
toroid is stretched toward the central object because of the strong 
gravity of the central object. In addition to this, the cusp-like shape 
at the inner edge of the toroids becomes 'sharper' as the value of 
m becomes smaller. We find the same tendency for the A'^ = 1.5 
polytropes. 

The top panel of Fig. |2] shows the structures of the magnetic 
fields on the meridional plane for the N — 3 critical configurations 
with m = —1.4 (solid curves) and m = 0.5 (dashed curves). In 
this figure, the surfaces of the toroids are indicated by the tear- 
shaped closure curves. The structures of the magnetic fields for the 
model with jn = — 1.4 are remarkably different from those for the 
model with m — 0.5. The shapes of the contours of the magnetic 
flux function for the toroid with m = 0.5 look nearly circle, but 
those for the toroid with m = —1.4 deform oblately. Moreover, 
the magnetic field lines are more densely distributed near the inner 
edge region of the toroid for the toroid with m = —1.4 compared 
to those for the toroid with m — 0.5. In other words, the magnetic 
fields are highly localized toward the central object for the toroid 
with m = —1.4 compared to those of the toroid with m — 0.5. 

The panels in the middle left and the bottom left of Fig. |2] 
show distributions of log^Q |B| on the equatorial plane, and the 
panels in the middle right and the bottom right of Fig.|2]show val- 
ues of each term in the right-hand side of the first line of Eq. (|7]l on 
the equatorial plane. Note that the horizontal axises of these figures 
range from Rinner/re{= qc) to Router/re{= 1). As seen from 
the middle left and the bottom left panels of Fig. |2] the distribu- 
tions of the magnetic fields are significantly different for the two 
equilibrium configurations. The ring of maximum magnetic field 
strength locates near the inner edge of the toroid for the model with 
m = —1.4, while it locates near the central region of the merid- 
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0.495 
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3.644E-05 


-1.1 


0.242 
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3.224E-05 


-1.4 


0.057 


1.467E-I-00 
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6.864E-01 


5.261E-01 


1.580E-01 


2.944E-05 



Table 1. Physical quantities for the critical configurations with Qo = 0.0 (no meridional flow), f!o = 0.0 (no rotation) and kq = 4.5. 
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Figure 1. The density contours on the meridional cross section for four critical configurations with different values of m. From left to right, the density 
contours correspond to the toroids with m = —1.4, m = —0.7, m = —0.1, and m = 0.5, respectively. Values of the other parameters are Af = 3, Qo = 0, 
Qq = 0, and kq = 4.5. The tear-shaped closure curves with a cusp-like structure indicate the surfaces of the toroids. The density difference between two 
adjacent contours is one-tenth of the maximum density. It is observed that the width of the toroid on the equator becomes wider as values of m decrease. 



ional cross section of the toroid for the model with m — 0.5. The 
toroids with m = —1.4 can sustain highly localized and strong 
magnetic fields in the nearer region from the central object com- 
pared to the toroids with m = 0.5 and can extend themselves 
toward the central object although the gravitational force of the 
central object is much stronger there. This implies that the toroids 
with m — —1.4 can produce strong magnetic force near their in- 
ner edge, which is balanced against the gravitational force of the 
central compact object as shown later. 

In the middle right and the bottom right panels of Fig. 



[2] the gravitational potential of the central object, — A/c/47rf 
(thick dashed curve), the gravitational pote ntial of the toroids, (j) „ 
(thin dashed curve), the magnetic potential I IFuiisawa et al.ll2012h . 
— J p. d'i/ — C (dotted curve) and the sum of all the potentials (thick 
solid curve) are shown as functions of Kjr^. Here, physical quanti- 
ties with (") are dimensionless quantities defined in Appendix |C2 1 
From these figures, it is clearly seen that the magnetic force is the 
primary agent supporting the toroid against the gravitational force 
of the central object. The gradient of the gravitational potential of 
the central object for the toroid with m — —1.4 is steeper than 
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Figure 2. Top panel: The poloidal magnetic fields on the meridional plane for the N = 3 toroids with m = —1.4 (sohd cui'ves) and rra = 0.5 (dashed 
curves). Values of the other model parameters are the same as those of the models given in Fig. [T] The tear-shaped closure curves with a cusp-like structure 
indicate the surfaces of the toroids. Middle left and bottom left panels: logj^Q |B| on the equator for two critical configurations with m = —1.4 (middle) 
and m = 0.5 (bottom), given as functions of R/re- The horizontal axises range from Rinner/'i'e{= Qc) to Routur / 'r e{= 1)- Middle right and bottom 

right panels: —Mc/i-K-r (thick dashed curve), (pg (thin dashed curve), — f fi('if) d'i — C (dotted curve), and sum of these three potentials (solid curve) on 



the equator for the critical configurations with m 

Rinner/re{= Qc) tO Router /re{= 1). 



-1.4 (middle) and m = 0.5 (bottom), given as functions of R/re- The horizontal axises range from 
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that for the toroid with m = 0.5 because the m — —1.4 toroid is 
located closer to the central object than the m — 0.5 toroid. The 
magnetic potentials behave very differently for these two equilib- 
rium configurations with different values of m. For the m — 0.5 
toroid, the magnetic potential curve has a substantial local mini- 
mum at R/ve ~ 0.9. For the m — —1.4 toroid, however, the 
magnetic potential curve is shallower and extends within a broader 
region and its slope is steepest near the inner edge of the toroid. 
As a result, the strong magnetic fields can exist near the inner edge 
region of the toroid and their magnetic force supports the toroid 
against the gravitational force of the central compact object. In this 
way, the m = —1.4 toroid can be in a stationary state even if the 
gravitational potential becomes much steeper as approaching to the 
central object. 



3.2 Effects of the meridional flows on the magnetized 
configurations 

3.2.1 Basic features of magnetized configurations with 
meridional flows 

Concerning the parameters which appear in the arbitrary functions, 
to examine effects of the meridional flow on toroid structures, the 
following values are chosen in this subsection: kq = 0.5, a — 
-0.5, d = 0.1, m = OandCll = 1.0 x 10"^ Note that smaller 
values of kq result in more rapid meridional flows and that this 
small value of Qq gives equilibrium configurations with almost no 
rotation. Parameters for the rotation law are the same as those in 
lYoshidaetai] ( l200d) . 

The left panel of Fig. [3] shows contours of the flux function VE" 
on the meridional plane for the critical configuration of an A'^ — 1.5 
polytrope with meridional flows. Direction of the fluid velocity on 
the meridional cross section in the critical configuration is shown 
in the right panel of Fig. [3] Here, the lengths of the vectors are not 
proportional to the absolute values of the fluid velocities. Note that 
the region where the meridional flows are present is only a part 
of the meridional cross section of the toroid because a particular 
functional form for Q'{9) is used to avoid singular behaviors of 
the meridional flow near the toroid surface (see Appendix IA21 . The 
bottom panel of Fig. [3] displays the velocity distributions normal- 
ized by the local Kepler velocity, on the equatorial plane. The solid 
and dashed curves denote the absolute value of the meridional ve- 
locity and the rotational velocity, respectively. As seen from this 
panel, the rotational velocity is sub-Keplarian, because our rota- 
tional parameter Qo is assumed to be small in this paper. On the 
other hand, the meridional velocity is slightly faster than the rota- 
tional velocity in this parameter region (Qo = 20). 

As shown in Figs. |2] and [3] (see, also. Table |2] given later), 
general structures of the toroids and their magnetic fields do not 
change significantly even when the fluid flows exist on the merid- 
ional plane. As argued later, however, the presence of the merid- 
ional flows changes the density distributions of the toroids slightly 
and increases the critical distance a little bit. 



3.2.2 Critical distances for magnetized toroids with meridional 
flows 

In order to study the influence of the meridional flows on the critical 
distances, we calculate two polytropic sequences (A*' = 1.5, A'^ — 
3) by changing the value of Qq for the m = toroidal current 
function. We fix the other parameters as kq ~ 0.5, a = —0.5, d — 
0.1, and Oo ~ 1-0 x 10^^ in this subsection. Physical quantities 



for several models belonging to the two polytropic sequences are 
tabulated in Table |2] Here, Tp denotes the kinetic energy of the 
meridional flow. 

As seen from this table, the energy ratio rp/|W| is much 
smaller than that of A1/|W^|, thus, the kinetic energy of the merid- 
ional flow is much smaller than the magnetic energy for the models 
obtained in the present study. This means that the magnetic forces 
mainly support the toroids against the gravitational forces of the 
central objects even when the meridional flow becomes stronger in 
the present parameter space. The meridional flow cannot change 
global structures of the toroids and their magnetic fields signifi- 
cantly. However, the energy ratio Tp/\W\ reaches ~ 1.0 x 10~^ 
when Qo > 40 (iV = 1.5) and Qo = 40 (W = 3). In fact, 
the critical distance Qc increases as Qo increases as shown in Ta- 
ble |2] This implies that the toroids tend to shed their mass when 
the meridional flow becomes stronger. In some sense, therefore, 
the influence of the meridional flow on the density distribution of 
the toroids should not be considered to be small. The rotational ve- 
locities of these models are also sub-Keplarian. The typical value 
is about 2 % of the Kepler velocity which is similar to that given 
in Fig[3] On the other hand, the meridional flow velocity is about 
several times as large as the rotational velocity. The maximum ve- 
locity of the meridional flow for the Qo = 80 model reaches about 
10 times as large as that of the rotational velocity. 

In order to clarify the effects of the meridional flows, let us in- 
vestigate the density distributions of toroids and the profiles of po- 
tential terms on the equatorial planes for the TV = 1.5 toroids with 
and without meridional flows. Fig. |4] shows the density contours 
on the meridional plane (top left and middle left) and the profiles 
of the potentials on the equator (top right and middle right). Bot- 
tom panels of Fig.|4]show the profiles of the density on the equator 
(left) and the velocity potential term, ||i)^| — Rvy^Cl, on the equa- 
tor (right) as functions of R/r^ . In the top right panel of Fig. |4] 
each curve denotes ^\v'^\ — Rv^tl (solid curve), ~Mc/4:-Kr (thick 
dashed curve), (j)g (thin dashed curve) and — J fid^ — C (doted 
curve). As seen from these profiles, in both the models, the gravita- 
tional potentials of the toroid make a tiny contribution to the equi- 
librium solutions and the balance between the teiTn — J /id\E' — C 
and the gravitational potential of the central object mainly deter- 
mines the stationary states of the toroid. Taking a detailed look at 
the middle right panel of Fig.|4] we observe that the potential terms 
due to the meridional flow and rotation (solid curve) show their 
maximum value near the radius of R/r^ = 0.95 for the A*' = 1.5 
models with Qo = 80.0. This very tiny protuberance is considered 
to appear due to the presence of the meridional flows because as 
shown in Table |2] the kinetic energy due to rotation is negligible 
small. More detailed structures of the velocity potential terms can 
be seen by enlarging these tiny protuberances. The bottom right 
panel of Fig.|4]shows the profiles of ||'&^| — Rv^tl on the equator 
for A'^ = 1.5 polytopes with Qo = 20 (solid curve), Qo = 40 
(dashed curve) and Qo = 80 (dotted curve). These profiles have 
double peaks which locate at i? ~ 0.85 and R ~ 0.92. These dou- 
ble peaks appear from the balance of the density distributions of the 
toroids and the meridional flows. As we have described in Sec. [T] 
the presence of the poloidal velocity fields results in reducing the 
density of toroids (see Eq. (|7}). For our numerical examples, the 
presence of the poloidal velocity fields decreases the density on the 
equatorial plane around radii of R ^ 0.85 and R ~ 0.92, which 
can be observed in the bottom two panels of Fig|4] (more detailed 
considerations are given below). 

The top left and middle left panels of Fig.|4]show the density 
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Figure 3. The poloidal magnetic fields (left), the meridional flow pattern (right) on the meridional plane for the toroid with N = 1.5, Qo = 20, Hq = 
1.0 X 10~^, m = 0, g = 0.6, and kq = 0.5. The tear-shaped closure curves with a cusp-like structure indicate the surfaces of the toroids. Length of 
vectors given in the right panel is not proportional to the speed of the meridional flows and does not have any physical meaning. The bottom panel shows the 
distribution of the absolute value of the meridional velocity (solid line) and the rotational velocity (dashed line) normalized by the local Kepler velocity on the 
equatorial plane. 



Qc 



An re (cm) m/\W\ 



n/\w\ 



T/\W\ 



Tv/\W\ 



N = 1.5 



N ■- 



0.625 2.810E-h00 
20 0.628 2.784E-H00 
40 0.632 2.711E-H00 



3.55 1E7 6.583E-01 
3.570E7 6.598E-01 
3.622E7 6.639E-01 



VC 



0.597 3.047E-H00 2.500E7 6.213E-01 1.262E-01 8.322E-05 O.OOOE-l-00 2.656E-05 

20 0.600 3.020E-H00 2.514E7 6.226E-01 1.255E-01 4.138E-04 3.185E-04 3.032E-05 

40 0.605 2.942E+00 2.555E7 6.268E-01 1.235E-01 1.431E-03 1.314E-03 4.271E-05 

60 0.615 2.817E-H00 2.627E7 6.338E-01 1.199E-01 3.217E-03 3.063E-03 6.505E-05 

80 0.629 2.653E-H00 2.759E7 6.437E-01 1.148E-01 6.101E-03 5.875E-03 1.595E-04 



1.133E-01 9.202E-04 O.OOOE-l-00 5.497E-05 
1.125E-01 1.464E-03 4.866E-04 3.642E-05 
l.lOlE-01 2.935E-03 1.842E-03 3.115E-05 



Table 2. Physical quantities for the critical configurations with meridional flows. Model parameters are kq = 0.5, fig = 1.0 X 10 ^ , and m = 0. 
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Figure 4. Top left and middle left panels: Density contours on the meridional cross section for the critical configurations with Qg = (top left) and Qq = 80 
(middle left). The heavy curves indicate the surfaces of the toroids. The density difference between two adjacent contours is one-tenth of the maximum density. 
Top right and middle right panels: ^\v | — Rv^Q. (solid curve), — Mc/47rf (thick dashed curve), — / A i^'I' — C (dotted curve) for the critical configurations 
with Qo = (top right) and Qo = 80 (middle right), given as functions of R/r^. The horizontal axis ranges from Rinner / "i" e{= Qc) to Router /Te{= !)■ 
Bottom left panel: Densities normalized by pc on the equatorial plane for the critical configurations with Qo = (dashed curve) and Qo = 80 (solid 
curve), given as functions of R/r^- The horizontal axises range from Rinner/fei= Qc) to Router/re{= 1). Bottom right panel: The dimensionless velocity 
potential term, ^ |i> | — Rv^Q, on the equatorial plane for the the critical configurations with Qo = 20 (solid curve), Qo = 40 (dashed curve) and Qq = 80 
(dotted curve), given as functions of R/r^. The horizontal axis ranges from Rinner/i"e{= <lc) to Router/re{= !)■ The model parameters are m = 0, 
N = 1.5, KQ = 0.5, and Qg = 1.0 X 10^^. 
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distributions of A'^ = 1.5 toroids with Qo = 0.0 and Qo = 80, re- 
spectively. Comparing the top left panel with the middle left panel, 
we see that the matter distribution of the toroid with Qo ~ 80 is 
shifted outward slightly. The bottom left panel of Fig|4]displays the 
density profiles on the equator for each toroid. The dotted and solid 
curves denote the density profiles for the models with Qo = 0.0 
and Qo = 80, respectively. As seen form this panel, the meridional 
flows shift the place where the density takes its maximum value 
outward and make the density gradient around R ~ 0.92 steeper. 
This is due to the double peak structure of the velocity potential 
profiles. The inner peak of this potential affects to decrease the den- 
sity around R ~ 0.85 where the density takes the maximum value 
if there is no meridional flow. On the other hand, the outer peak of 
this potential also leads to decrease in the density around R ~ 0.92. 
Since the density decreases around i? ~ 0.85 and _R ~ 0.92 by the 
presence of the meridional flows, the place where the density be- 
comes maximum moves outward and the density gradient becomes 
steeper if the toroids have rapid meridional flows (Qo = 80 model). 
These effects result in decreasing the critical distances. 

Next, we deal with the influence of the equation of state on 
the critical distance. As we have seen in Table |2] the critical dis- 
tances of the N — 3 toroids are larger than those of the A'^ = 1.5 
toroids. The same tendency in the equilibrium c onfigurations with- 
out meridional flows found in lotani et al.l ( 12009.) . This is because 
that the mass shedding from the inner edge of the toroids is more 
likely to occur for softer equations of state. 

3.2.3 Effects of meridional flows on equilibrium conflgurations 
with highly localized poloidal magnetic flelds 

Finally, we unveil the effects of the meridional flows on structures 
of the toroid having highly localized poloidal magnetic fields. We 
consider the A'^ = 1.5 toroid models only in this subsection be- 
cause basic properties are independent of the equation of state. Fig. 
|5] displays typical models characterized by highly localized mag- 
netic fields with and without strong meridional flows. Here, we take 
m = —1.4, with which highly localized poloidal magnetic fields 
are obtained inside the toroid as argued in Sec. |3] Other parame- 
ters are kq = 0.5, a = -0.5, d = 0.1, and Clf-, = 1.0 x 10"•^ 
which are the same as those used in Sec. 13.2.2] The top two panels 
of Fig.[5]show the density distributions on the meridional cross sec- 
tion for the toroids with no meridional flow (left) and with strong 
meridional flows (right). The bottom two panels of Fig. |5] display 
the densities (left) and the velocity potential terms (right) on the 
equatorial plane as functions of R/re- In the bottom left (right) 
panel, the solid and dashed curves are correspond to the models 
with Qo = 15 (Qo = 15) and Qo = (Qo = 10), respectively. 
Comparing the bottom two panels of Fig. |5]to those of Fig.|4] we 
observe that the density profiles of the jn = — 1.4 models are sub- 
stantially different from those of the m = Q models though behav- 
iors of their velocity potential terms are similar in the sense that 
they show similar double peaks. The toroids with highly localized 
magnetic fields (the m — —1.4 models) are extended inward due 
to the strong gravity of the central object in comparison with the 
m — models (Compare the bottom left panels of Figs. |4] and |5] 
). As shown in the bottom left panel of Fig.|5] the positions of the 
maximum density rings for the m = —1.4 models are nearly inde- 
pendent of values of Qo, but the density gradient of the Qo — 15 
model is steeper than that of the Qo — model around R ~ 0.85 
because of the presence of the meridional flows. As a result, the 
matter distribution around _R ~ 0.7 is stretched outward. The pres- 
ence of the meridional flows also decreases the critical distance 



slightly as shown in the bottom left panel of Fig. \5\ which can be 
seen more clearly for the m = models (see the bottom left panel 
of Fig.©. 

As we have exhibited in numerical examples so far, values of 
Qc decrease as values of m decrease, while they increase as val- 
ues of Qo(> 0) increase. In other words, the poloidal magnetic 
fields generated by positive toroidal cuiTents are apt to expand the 
toroids to the directions normal to the equi-flux function surfaces in 
particular when the magnetic fields are highly localized around the 
inner edge of the toroids and the meridional flows act as an agent 
for shrinking the region where the fluid matter occupies. In order to 
quantify the influence of the highly localized magnetic fields and 
the meridional flows on the critical distance, we introduce a quan- 
tity Ag^(Qi, mi; Q2,m2), defined by 



A<j^(Qi,mi;Q2,m2) = 



qc{Qi,mi) 
qc{Q2,m2) 



(24) 



where qc{Qo,tno) is the critical distance of the equilibrium se- 
quence of the toroid characterized by A'^ = 1.5, Alt/Mc ~ 
2.0 X 10~^, Q — Qo, and m — mo. Positive (Negative) val- 
ues of A,^ mean that the critical distance for the sequence with 
(Qi, mi) increases (decreases) or its width of the toroid decreases 
(increases) compared to that with (Q2, m,2). In the left and right 
panels of the Fig. |6] A,^ (0, m; 0, 0) is given as a function of m 
and Aq^{Qo,Tn; 0.,m,)'s are given as functions of Tp/jVt^j for sev- 
eral fixed values of m, respectively. The left panel shows that the 
values of Ag^ range from about —0.9 to 0.2. Highly localized mag- 
netic fields (for models with negative values of m) show the signifi- 
cant influence on the critical distances. On the other hand, the right 
panel shows that A,^ can reach about 0.05 due to the effects of 
the meridional flows. Regardless of the sign of m, the values of 
Tp/I W^l range up to ~ 0.006 and the maximum value of A,^ can 
reach 0.05 as the values of Qo are changed. Thus, the maximum 
value of A,^ would be 0.05 when the meridional flows exist. This 
means that the influence of the meridional flows on the critical dis- 
tances is much smaller than that of the highly localized magnetic 
fields, but it is certainly true that the meridional flows work as an 
increasing factor for the critical distances. We also find that the ef- 
fects of the poloidal magnetic fields and the meridional flows may 
nearly cancel out for the toroids characterized by Tp/j PFI ~ 0.005 
and m = —0.2. For this model, the critical distance or the width 
of the toroids is similar to that of the model with Tp/\W\ = 
and m = 0. As expected in Sec. 11.11 thus, we confirm that the 
oppositely working effects of the highly localized magnetic fields 
and the meridional flows result in nearly no change in the critical 
distance for some particular toroid model. 



4 DISCUSSION AND CONCLUDING REMARKS 

4.1 Strength of the magnetic flelds inside the toroids 

As shown in the middle left and bottom left panels of Fig. [2] for 
some particular set of the parameters, the toroids in the critical 
states have very strong magnetic fields and their strength are about 
lO^'^ G if we take Mc = 5.0Mg and pc = 1.0 x 10"g/cm^ 
Fig. |7] shows the magnetic energy of the toroids as functions of m. 
Here, the magnetic energy is normalized by that of the m = 0.0 
and A' = 3 model. From this figure, it is found that the magnetic 
energy becomes larger as the value of m decreases. The magnetic 
energies of the m = — 1.4 models are nearly three times larger than 
that of the m = 0.0 and A'^ = 3 model. Therefore, we conclude that 
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Figure 5. Top panels: The density contours on the meridional cross section for the critical configurations with Qq = 0.0 (left panel) and Qg = 15 (right 
panel). The heavy curves denote the surfaces of the toroids. The density difference between two adjacent contours is one-tenth of the maximum density. Bottom 
left panel: Densities normalized by pc on the equatorial plane for critical configurations with Qo = (dashed curve) and Qo = 15 (solid curve), given as 
functions of R/re- The horizontal axis ranges from Rinner/i"e{= Qc) to Router/re{= 1). Bottom right panel : The dimensionless velocity potential term, 
i I'D I — Rvipfl, on the equatorial plane for critical configurations with Qo = 10 (dashed curve), Qo = 15 (solid curve), given as functions of R/re- The 
horizontal axis ranges from Rinner/i'e{= Qc) to Router/'>'e{= 1). The model parameters are m = —1.4, A'^ = 1.5, kq = 0.5, and Q^ = 1.0 X 10~^. 
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Figure 6. Aq^(0, m, ;0, 0) vs. m (left panel) and Aq^(Qo,"^i ; 0, m) vs.Tp/|Vy| for the critical configurations with m = 0.5, 0, —0.7, —1.4 (right panel). 
Here, Ag^{Qi, mi; Q2,"i2) = qdQi, 'nT-l)/<lc(Q2, "^2) — 1 and Tp/ jiyj is given as a function of Qo for the sequence of the critical configurations. 



larger magnetic energy can be sustained in the toroids whose mag- 
netic fields are highly localized around the inner edge of the toroids 
and which locate closer to the central compact object. 



4.2 Critical distances for magnetized toroids with meridional 
flows and highly localized magnetic fields 

Assuming the toro idal current function fi to be constant (m — 0), 
lOtani et alj ( 120090 showed that there appears a critical distance in 
the self-gravitating toroids with the magnetic fields and that the 
critical distances are much larger than the radii of the inner-most 
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Figure 7. The magnetic energy normalized by the value of the equilibrium configuration with m = 0.0 and N = 3, given as functions of m. The solid and 
dashed curves correspond to the A'^ = 1.5 and N = 3 polytropes, respectively. 



Stable circular orbit (ISCO) of the Schwarzschild black hole with 
the mass M = Ale- In this study, we show that their conclusions 
hold true when the meridional flows are taken into account and the 
functional form of /i is generalized to the cases of 7n > 0. For the 
cases of 771 < 0, however, we find that the critical distance can be 
much shorter than that of the m = case. The radius of ISCO of 
the Schwarzschild black hole with the mass AI = Mc, risco, is 
given by 



nsco 



6GMc 



(25) 



If Adc — 5.0Mq, which is the fiducial value in this study, risco ~ 
4.43 X 10® cm. For the A'' = 3 models with m = 0.0 and m = 
— 1.4, using results given in Table [T] we respectively obtain 



nn = g X re 



2.30 X lO^cm > nsco : 



and 



g X re ~ 1.08 X 10 cm < nsco ■ 



(26) 



(27) 



Therefore, the critical distance can be the same order or even 
smaller than the radius of the ISCO in this parameter space. Since 
this study is done within the framework of Newtonian gravity, the 
quantitative evaluation is not correct if rin — O(risco), while 
the results given in the present study are reasonable as long as 
rin 2> risco- We need to use general relativity for the toroids with 
rin ~ O (risco)- However, this is beyond the scope of the present 
study as mentioned before. 

4.3 Concluding remarlss 

In this paper, we have investigated and calculated the stationary 
states of magnetized self-gravitating toroids with meridional flows 
and various kinds of the magnetic fields around central compact ob- 
jects. As a result, we have obtained the toroids with strong merid- 
ional flows and with strong poloidal magnetic fields. Our findings 
and conjectures are summarized as follows. 

(i) Choosing the functional forms of arbitrary functions, we can 
change the strengths of the meridional flows. The critical distances 
for stationary toroids with meridional flows become larger than 
those for stationary toroids without meridional flows. In addition 
to this, the distances increase as the strengths of meridional flows 



becomes larger. This is what we have discussed in Sec.[T] i.e. the ef- 
fects of the meridional flows and the magnetic fields are oppositely 
working. 

(ii) Changing the value of the parameter m in the certain choice 
of the arbitrary function /i, we can change the distributions of the 
poloidal magnetic fields inside the toroids. In particular, the crit- 
ical distances could be smaller as the value of m is decreased. If 
we adopt Ale ~ 5.0Mq and the maximum density of the toroid 
to be /Omax ~ 1-0 X lO^'^g/cm , the critical distance for the 
m = —1.4 toroids becomes the same order as that of the ISCO 
of the Schwarzschild black hole with mass M^. For such toroids, a 
general relativistic treatment is necessary for their correct descrip- 
tion. 

(iii) The magnetic energy for the critical configuration could in- 
crease as the value of m is decreased. The magnetic energy for the 
critical configuration with m — —1.4 is about three times larger 
than that for the m — 0.5 critical configuration. 

(iv) By obtaining stationary configurations of axisymmetric 
magnetized self-gravitating polytropic toroid with meridional flows 
under the ideal MHD approximation, we have shown that the ef- 
fects of the meridional flows would work oppositely to those of 
the poloidal magnetic fields. In other words, the oppositely work- 
ing effects can be easily understood if we consider that the dense 
magnetic field lines expand the gaseous configurations due to the 
repulsive nature of the magnetic field lines and that the presence of 
the meridional flows works as lowering the gas pressure due to the 
appearance of the ram pressure as seen from the stationary condi- 
tion equation. 
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APPENDIX A: INTEGRABILITY OF BASIC EQUATIONS 
AND APPEARANCE OF ARBITRARY FUNCTIONS: 
MAGNETIC FLUX FUNCTION BASED FORMULATION 

A I Integrablllty conditions and expression for the current 
density 

Since equilibrium configurations treated in our problem are ax- 
isymmetric, Eqs.([8} and ill) are automatically satisfied by intro- 
ducing the stream function Q{R, z) and the magnetic flux function 
*(7?,2), defined by 

I dQ _ 1 dQ 



VR 



and 



47rp_R dz 



A-KpR dR 



(Al) 



1 a* „ 19* 
^- = -^97' ""^^RdR- ^^2^ 

If and only if the meridional flows and the poloidal magnetic fields 
exist simultaneously, the stream function is given by a function of 
the magnetic flux function, i.e., 

Q = Q(*). (A3) 

This relation is obtained from Eqs. dlS) and ( |15t . From the ip- 
components of Eqs.lfTsll and l ll5t . we obtain 



~A ^f 

A-Kp 



R 



= ^W , 



(A4) 



where means the derivative with respect to '^ and f2(*I') is an 
arbitrary function of "if. From integrability conditions for Eq. (|9), 
we obtain two relations: 



RB^ - Q'(*)i?u^ = k(*) , 



(A5) 



AwpR 



Rv^n'{<i!) ■ 



Vv>Q"i'^)B^ 
Aizp 



pR 



?'(*) 



4wpR 



M*). (A6) 



where k(*) and p{'i') are arbitrary functions of *. Using Eqs. 
l ll4t and l|A6), we may describe the current density j as follows: 



I 

c 



+pR [At(*) + Rv^n'i'i')] e^ . 



(A7) 



It should be noted that in this current density formula, the four 
arbitrary functions appear and that each function is related to differ- 
ent physical qua ntities corresponding t o each term in the right-hand 
side of Eq. ( IA7t . lLovelace et al.l ( ll986f) obtained the current density 
including another arbitrary function which is related to the entropy. 
This arbitrary function, related to the entropy, does not appear in 
our problem because uniform entropy distributions is implicitly as- 
sumed. 



A2 Choice of arbitrary functions in this study 

We choose the functions k, Q' , fl and p as: 

, for * < *niax , 



k + 1 



(* - *max)* + ' , for * > *„ 
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Q'(*) = 
n(*) = 



fc + i 



(*-eQ*n 



U+1 



for 
for 



* < eg*!, 



where the self energy of the central object is discarded. The kinetic 
\"^) energy of the fluid is defined by 






(AlO) 
(All) 



where, ko, k, Qo, eg, ^o, d, fio, e, a, and m are con- 
stant parameters and ^'^ax denotes the maximum value of ^ 
in the vacuu m region. T his choice of k i s the same as that 
emplo y ed in Tomimura & Eriguchil ( 20051). Yoshida & Eriguchil 
( l2006h . lYoshida et alj ( l2006h . and lOtani et all J2009l) . We introduce 
the parameter eg in Eq. JA9b and set eg = 1.1 in order to restrict 
the region where meridional flows exist well inside the surface of 
the toroid. Choosing these functional forms, we can avoid singular 
behavior of the solutions which could appear on the surface of the 
toroid. 



APPENDIX B: STREAM FUNCTION BASED 
FORMULATION: VORTICITY FORMULA 

So far, we consider the situations in which the poloidal magnetic 
fields exist everywhere except in vacuum region. For such situ- 
ations, as already shown, the flux function ^ can be a principal 
variable by which all the magnetohydrodynamical quantities are 
determined. If we assume that the poloidal velocity fields exist ev- 
erywhere inside the fluid region, the same problem as that treated in 
this study may be formulated by considering the stream function Q 
as a principal variable, which is named "the stream function based 
formulation". For this formulation, the magnetic flux function is 
given by a function of the stream function. 



* = *(Q) 



(Bl) 



The other arbitrary functions of the magnetic flux function defined 
in Appendix lAl are regarded as functions of Q. Then, the vorticity 
vector may be written as: 



47rp 



d£{Q) ^ ^ d^^jQ) ^^ 



+ 47r 



dQ 
dQ 



dQ^ 
'-+pR\-v{Q) 



jjdaiQ) 



(B2) 



where £{Q), v{Q) and cr(Q) are another set of arbitrary functions 
of Q for the stream function based formulation. The arbitrary func- 
tions £{Q) and (t(Q) are related to the physical quantities Vip and 
By, through 



£{Q) = Rv^ 






B^ 47r d^(Q) 



(B3) 



(B4) 



2 



p\v\ d r . 



The volume integral of the pressure is defined by 

n= pd^r . 

Through H, the internal energy is defined as 

U = NU. 

The magnetic energy is defined by 



OR; 



- X B 



j3 

a r . 



The mass of the toroid is defined by 



M 



pd r . 



(C2) 



(C3) 



(C4) 



(C5) 



(C6) 



In order to evaluate the effect of the meridional flows, we also de- 
fine the kinetic energy of the meridional flows as follows: 



Tr,= 



pVp 



2 ,3 

a r. 



(C7) 



C2 Dimensionless physical quantities 

For the numerical computations, the physical quantities are trans- 
formed into dimensionless ones. The dimensionless quantities em- 
ployed in this study are defined as follows: 
For local quantities. 



Te 



V ;9 47T 



G^ 



^ 4% Gripe 



fi; 



V47rGpc 



^ATvGr'ipc ' 



viTrGrlpc 
^4^/re ' 



(C8) 

(C9) 
(CIO) 
(Cll) 
(CI 2) 
(C13) 
(C14) 



APPENDIX C: PHYSICAL QUANTITIES 

CI Global physical quantities 

Some useful global quantities are defined as follows: The gravi- 
tational energy of the magnetized toroid and the central object is 
defined by 



^=2 



GA/c 



pd'^r , 



(CI) 



\/47rGrePc 



*; 



v47rGr|pc 



V47rGr|pc 



(C15) 



(CI 6) 



(C17) 
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K : 



J-fi 



c 



K 






V47rGpcC 
C 



(CI 8) 



(C19) 



(C20) 



47rGripc ' ' 

where r^ and Pmax denote the equatorial radius of the outer edge 
of the toroid and the maximum pressure Pmax, respectively. 
As for the global quantities, 

M 



M 



W 



rlp^ 



W 



A,-KGr^,pi ' 



(C21) 



(C22) 



vacuum 




► R 



R=q / R= 1 



P = Pc = i 

Figure Dl. Three points used in the HSCF scheme. 



n = 



m 



n 

4nGr^pc 



(C23) 



(C24) 



(C25) 



47rGrip2 ■ 

Here, /3 appearing in Eq. dCSt is introduced so as to make the dis- 
tance from the symmetry axis to the outer edge of the toroid to be 
unity (or f = 1) during numerical iterations (see, e.g.. lOtariiet al.l 
1 120090 ). Thus, values of /? are obtained after converged solutions 
are obtained. Increasing /3 results in decrease in r^ and decreasing 
(3 results in increase in r^- 

Using the dimensionless quantities defined above, we obtain 
the dimensionless form of the right-hand side of Eq. (|7j. 



AnGr'ipc J p 47rf 2 

+ f sin 9v^Cl{i') + C 






(C26) 



We integrate the left side of this equation by using polytrope rela- 
tion as follows: 



47rGr|pc 



dp 



■P{N + I)p^^''. 



(C27) 



Then, we obtain the dimensionless form of Bernoulli's equation as 
follows: 



I3{N + I)p 



-l/N 



+ fsin6l{}^0(*) + C. (C28) 



We have used this dimensionless forms in actual numerical compu- 
tations. 



APPENDIX D: NUMERICAL METHOD 

In our numerical studies, the generalized iteration scheme known as 
Hachisu's Self-Consistent-Field (HSCF) method (Hachisu 1986) is 
adopted in order to solve the system of non-linear partial differen- 
tial equations for equilibrium configurations of magnetized toroids. 
In this generalized HSCF method, the density, the gravitational po- 
tential and the vector potential are discretized on grid points (r; , 9j ) 



and non-linear algebraic equati ons for these discretized ph ysical 
quantities are iteratively solved (ITomimura & EriguchilBOOSl) . 

We start our computation by assuming initial guesses for the 
discretized mass density and flux function. Using a trial density 
and a flux function distribution, the dimensionless gravitational po- 
tential <f)cj and the vector potential A^p are respectively computed 
through the integrations, which are equivalent to Eqs. ( 120 ) and J21l >. 



h{f-,d) = -^P2n(cOSe)/ df'f'^f2-a{r,f') 
n-K/2 

X / de'sin6''P2„(cos6l')p(r',6l'), (Dl) 

Jo 



A,{r,e) = ^^''f^^lr.-.{cosi 



^ 2n{2n - 1) J^ 



df'r /2n-i(r, f') 



X / de'sme'Pl„_^{cose')j^{r',9'), (D2) 
Jo 



f'2n{f,f') 



~2n /^/2n + l 



{f>f'), 
(f < f'). 



(D3) 



Here P„(cos&) are the Legendre polynomials of degree n and 
Pji(cos6') are the associated Legendre functions of order one. In 
the actual computations, we take rimax = 40 typically. We use Eq. 
l |A7t to obtain the distribution of j^p. After these integrals, we cal- 
culate the distribution of the flux function using the relation of Eq. 
l ll6t . Although there are many constant parameters in our formula- 
tion, three constants, the length scale factor /3, the integration con- 
stant G, and the strength of the toroidal current po are especial in 
the sense that they cannot be specified before converged solutions 
are obtained. These three constants are determined by imposing 
the following three conditions given at three special grid points and 
therefore change their values during the iteration procedures. Fig. 
IDll shows a schematic image of our scheme. At the inner surface 
(point Q), i.e. at f = g and 9 = 7r/2, and at the outer surface (point 
P), i.e. at f = 1 and 6 — -k 12, the density must vanish. At the 
point where the density takes its maximum value (point C), i.e. at a 
point of f = Tmax and Q = ^max, the dimensionless density must 
be unity. From Bernoulli's equation I IC281 and three conditions are 
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explicitly given by 
= 



P + 



47rf I 



;{€■ +ye +^5^)1^ 



Mo 
m+ 1 



(*|p + e)™+^+fsin6»fi^no(* 



-d')"ip +(7,(04) 



n ill -'^'^ I -'- /"2 I -2 I -2m 



+ 



MO 
m + 1 



(*Iq + ^) 



m + l 



+ f sinei)^no(*^ + d^riQ + C , (D5) 



A4c 



/?(l + iV) = -0,ic + -4L 
47rr Ic 



-{vr +V0 +v^)\c 



+ ■ 



MO 



Wc + ey"+' + fsmOv^Qoi'i' + d^ric + C , (D6) 



Solving these equations, we obtain three constants /3, C and /io- 
Using the three constants, gravitational potential and flux function 
obtained newly, we solve Bernoulli's equation in terms of the mat- 
ter density as follows: 



(no rotation) for several different values of Nr- Fig. lEll shows val- 
ues of log VC against values of log Nr . In this figure, we observe 
that as the grid number in the r-direction is increased, the value of 
VC becomes smaller as 



VC oc ivr 



(E3) 



This behavior of VC is quite reasonable and proper because we use 
the numerical scheme with the second-order accuracy formula for 
numerical derivatives and Simpson's integration formula. If Nr = 
513 (log Nr ~ 2.710), the value of VC is ~ 10~^ (log VC = 
—4.58). In the present study, we choose the grid number as Nr = 
513 and No = 513 because this grid number is sufficiently large 
to obtain numerical solutions with acceptable accuracy as shown 
before. 



This paper has been typeset from a TgX/ ETgX file prepared by the 
author 



47rf 



{v^+vj + v'^) + ;|5^(* + e)'"+i +rsin0i)^Qo("I'^ + d^)°' + C ^N 



llf,-i 



P(l + N) 



(D7) 



The newly obtained density and other quantities are used as 
a new guess for the next potentials in the iteration cycle. We carry 
out this iteration procedure until the relative changes of all physi- 
cal quantities between two iteration cycles becomes less than some 
prescribed small number, e.g., 10^ . 



APPENDIX E: NUMERICAL ACCURACY CHECK 

In order to check the convergence and the accuracy of the obtained 
stationary solutions, we use a virial relation, 



2T + W + ?,n + m^Q. 



(El) 



This relation is satisfied for the exact stationary solutions. Numer- 
ical solutions, however, have numerical errors. Thus, in general, 
the right-hand side of Eq. dElt does not vanish although its value 
should be sufficiently small. In order to quantify numerical errors 
contained in the numerical solutions obtained, we evaluate a virial 
relation, defined by 



VC = 



2T + W + 3n + OT 



W 



(E2) 



VC is an estimate o f th e numerical errors and has been used 
in many papers (e.g.,lHac hisu 1986; Tomimura & Eriguchi 20051: 
iLander & Jonesll2009l : lotani et al.ll2009l:lFujisawa et al.ll2012h . If a 
value of VC is less than lO""* - 10~*, the solution is considered 
to be sufficiently accurate. We have checked the convergence of 
VC by changing the gird number in the f direction, Nr, to test 
our numerical code in which we have used the spherical coordi- 
nates. As for the grid number in the 9 direction, A''^, we fix it as 
No — 129 during the test computations. Since we assume equa- 
torial symmetry and the gravitational potential of the central ob- 
ject has a singularity at f = 0, we consider the region defined 
by f G [1.0 X 10"^ : 4.0] and 9 £ [0 : 7r/2] as our compu- 
tational space. In order to resolve structures of the t oroid properly, 
we us e non-uniformly distributed grid points used in lFuiisawa et al.l 
l2012h . We compute equilibrium configurations with q = 0.6, 
Mt/Mc = 2.0 X 10"^ m = 0,ko= 4.5, Qo = 0.0, fio = 0.0 
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Figure El. Values of VC vs. the number of the grid points in the r direction. 



